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ABSTRACT 

In the next years, several cosmological surveys will rely on imaging data to estimate the redshift of 
galaxies, using traditional filter systems with 4—5 optical broad bands; narrower filters improve the 
spectral resolution, but strongly reduce the total system throughput. We explore how photometric 
redshift performance depends on the number of filters n/, characterizing the survey depth by the 
fraction of galaxies with unambiguous redshift estimates. For a combination of total exposure time 
and telescope imaging area of 270 hrs m 2 , 4—5 filter systems perform significantly worse, both in 
completeness depth and precision, than systems with nt > 8 filters. Our results suggest that for low 
rif, the color- redshift degeneracies overwhelm the improvements in photometric depth, and that even 
at higher n/, the effective ■photometric redshift depth decreases much more slowly with filter width than 
naively expected from the reduction in S/N . Adding near-IR observations improves the performance 
of low n/ systems, but still the system which maximizes the photometric redshift completeness is 
formed by 9 filters with logarithmically increasing bandwidth (constant resolution) and half-band 
overlap, reaching ~ 0.7 mag deeper, with 10% better redshift precision, than 4 — 5 filter systems. A 
system with 20 constant- width, non-overlapping filters reaches only ~ 0.1 mag shallower than 4 — 5 
filter systems, but has a precision almost 3 times better, Sz — 0.014(1 + z) vs Sz — 0.042(1 + z). We 
briefly discuss a practical implementation of such a photometric system: the ALHAMBRA survey. 

Subject headings: Cosmology; photometric redshifts; galaxy surveys 



1. INTRODUCTION 

Photometric redshift estimation is not a new technique 
(Baum 1962, Loh & Spillar 1986, Koo 1985, Connolly et 
al. 1995; see Koo 1999 for a history of the method), 
but it has considerably developed in the last decade, es- 
pecially following the Hubble Deep Field (HDF) obser- 
vations (Williams 1996, Casertano et al. 2000), which 
provided catalogs with excellent photometric quality and 
abundant spectroscopic redshift coverage. This allowed 
astronomers to thoroughly test standard photo-z tech- 
niques and try new approaches (Gwyn & Hartwick 1996, 
Lanzetta, Yahil & Fernandez-Soto 1997, Sawicki, Lin & 
Yee 1997, Fernandez-Soto, Lanzetta k, Yahil 1999, Brun- 
ner et al. 1997, Beni'tez et al. 1999, Bem'tez 2000, Bol- 
zonella, Miralles & Pello 2000). 

As Hickson, Gibson & Callaghan (1994) first showed, 
multiband narrow filters can be much more efficient to 
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obtain redshifts than spectroscopy if the large area of the 
imaging cameras is factored in. Several photometric sur- 
veys, using different filter systems, have been proposed 
or implemented in the last decade: the UBC-NASA sur- 
vey (Hickson & Mulrooney 1998), CADIS (Wolf et al. 
2001b), COMBO-17 (Wolf et al. 2001a), COSMOS-21 
(Taniguchi 2004), ALHAMBRA (Moles et al. 2008), 
DES (DES collaboration 2005), LSST (Tyson 2006), 
PanStarrs (Kaiser 2007), VST (Arnaboldi et al. 2007), 
and PAU (Bem'tez et al. 2008). These surveys repre- 
sent powerful alternatives to deep spectroscopic surveys 
like DEEP2 (Davis et al. 2003), VVDS (Le Fevre et 
al. 2003), or BOSS (Schlegel et al. 2007) at least for 
those scientific goals which only require limited redshift 
accuracy and low resolution spectral information. 

However at least three of the imaging surveys (DES, 
LSST, PanStarrs) will work with photometric systems 
with 4 — 5 optical broadband filters, similar to those tra- 
ditionally used in Astronomy It is obvious that using 
more, narrower filters, for a fixed exposure time, will 
significantly sacrifice photometric depth. However, pho- 
tometric depth is not equivalent to photometric redshift 
depth. The fewer the filters, the more prone the system is 
to color-redshift degeneracies; these make impossible to 
unambiguously determine the redshift for a galaxy, even 
if observed at relatively high S/N. The Hubble Ultra 
Deep Field (Beckwith et al. 2003) offers a good exam- 
ple. Despite the fact that the limiting magnitude in the 
HUDF is 0.9 - 1.4 mag deeper than the HDF, the lack 
of a U— band filter in the HUDF makes the photometric 
redshift depth of both fields similar (Coe et al. 2006). 

This letters explores the impact on photometric red- 
shift performance of factors as the number of filters n / , 
constant vs. logarithmically increasing bandwidth, half- 
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Fig. 1. — Example of the four types of filter sets considered, each 
with 11 filters. We represent the filter transmissions without taking 
into account the CCD or the telescope+optics transmission, factors 
which are later included in the photometric noise estimation to 
produce realistic photometric measurements. We slightly increase 
the height of the filters with wavelength, and alternate colors in 
successive filters to help visualization. 

Examples of filter systems 
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band overlaps, and near-IR observations. We also briefly 
discuss a practical implementation of a medium band fil- 
ter system: the ALHAMBRA survey. 

2. SIMULATIONS 

2.1. Description of the filter systems 

We assume that filters are almost "top-hat", with a 
transmission which is constant in the central part and 
steeply falls on the edges, formed by half-Gaussian wings 
with a HWHM of ~ 15 A. Although somewhat ideal- 
ized, this is very similar to the characteristics of the fil- 
ters provided by BARR Associates for the ALHAMBRA 
survey. We consider four type of photometric systems, 
depending on whether they have constant or logarithmi- 
cally increasing (AA ~ A) bandwidth, and whether they 
have half- width AA/2 overlaps or just a minimal overlap 
corresponding to the filter wings. The filters cover the 
3400 — 9800Ainterval. Fig. 1. shows examples of the 4 
types of filter systems considered. 

2.2. Mock catalogs 

To generate realistic galaxy distributions, we use the 
HDF catalogs of Fernandez-Soto, Lanzetta & Yahil 
(1999) and Yahata et al. (2000). It is possible 
to obtain accurate and reliable photometric rcdshifts 
and type classifications, based on the extended Cole- 
man, Wu and Weedman (1980) set for these galax- 
ies up to Iab ~ 27 (Bem'tez 2000, Fernandez-Soto et 
al. 2001). We use the distribution of Iab, Bayesian 
spectral type if, and Bayesian photometric redshift Zb 
obtained from these catalogs with the software BPZ 
(Ihttp: / / acs.pha.jhu.edu/~txitxo , Bem'tez 2000) and the 
empirically calibrated spectral library of Benftez et al. 
2004. The original input catalogs contain 822 galaxies 
with Iab < 26. For each filter set combination, our sim- 
ulation includes 5,000 galaxies, generated by randomly 
choosing objects from the HDF catalogs. The depth 
reached by our simulations (up to / ~ 25) precludes us- 



ing an input catalog based on spectroscopic redshifts, 
which would not be complete enough at those magni- 
tudes. Since the accuracy of the input photometric red- 
shifts is w 0.06(1 + z), we perturb them by a similar, 
randomly distributed amount to produce a more realis- 
tic redshift distribution. 

We assume a constant total exposure time of 100, 000s 
(28hrs) per pointing, and vary the number of filters rif. 
The average exposure is thus 100,000/n/. Instead of 
using a fixed observing time per filter we distribute the 
exposures trying to reach constant S/N for a same AB 
magnitude, with two constraints: the minimal exposure 
time per filter is, for practical reasons, at least 2, 500s, 
and we do not expose more than twice the length of the 
average exposure, i.e. 2 x 100,000/n/ to avoid spending 
too much of the total time on the filters which are less 
efficient. For wavelengths A > 8000Athe efficiency is so 
low that we set this limit to be 100,000/n/. 

Using this information, we generate the expected mag- 
nitudes in all the filter systems described in the above 
section using functions included in the BPZ package. For 
~ 300A filters the corresponding 5a limiting magnitudes 
are ttiab ~ 25 bluewards of 8000A and quickly degrades 
to rriAB ~ 23 at 9500A. As it is obvious, for a fixed 
total exposure time, the limiting flux in a filter will be 
roughly inversely proportional to the square root of the 
filter width A A. 

As a last step we add random photometric noise whose 
amplitude has been estimated using the WHT exposure 
time calculator 11 scaling it to a 3.5m telescope. The 
product of exposure time by telescope area is ~ 270 hrs 
m , and therefore the results obtained here are equivalent 
to e.g. what could be obtained in a few hours of total 
exposure time with a 8 — 10m telescope. 

We also scale the S/N as a function of the magnitude of 
the galaxies. We do not take into account the dependence 
of size, etc. with magnitude. These are second order 
effects which will have a similar impact on all the filter 
systems and therefore are not expected to significantly 
affect the comparison among them. 

3. COMPARISON AMONG DIFFERENT FILTER SYSTEMS 

The photometric redshifts for the mock catalogs are es- 
timated using the BPZ package. The software provides 
a Bayesian estimate of the redshift and a spectral type 
classification. The expected reliability of the photo-z can 
be gauged through the Bayesian odds. The value of this 
parameter corresponds to the amount of redshift proba- 
bility concentrated on a ±0.2(1 + z) region around the 
probability maximum. Low values of the odds indicate a 
multimodal or very extended, little informative p(z), in- 
dicating that the photometric information is insufficient 
to obtain a unambiguous estimate of the galaxy redshift. 
By selecting objects with high odds, e.g. > 0.99, one 
can produce highly reliable samples (Benftez 2000), with 
very good redshift accuracy and a very low rate (< 2%) 
of "catastrophic" outliers. Therefore, it is possible to ac- 
curately characterize the effective completeness of a pho- 
tometric redshift catalog by using the amount of galaxies 
with odds above a certain threshold, which tells us how 
many galaxies we can expect to have meaningful, univo- 
cal photometric redshifts. 

11 Based on the SIGNAL software, by Chris Benn 



Optimal filter systems for photo-z 
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Fig. 2. — Effective 80% completeness magnitude, corresponding 
to the magnitude at which the accumulated number of objects N(< 
"10.99) with Bayesian odds > 0.99 is 80% of the total number of 
objects N(< m), a good measure of the effective depth of a survey. 
The blue dotted line illustrates how the completeness magnitude 
would change with filter number if it mimicked the behavior of the 
photometric limiting magnitude. 
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For a set-up with a total exposure time T and to- 
tal number of filters rif, the signal-to-noise in an indi- 
vidual filter, assuming that we are limited by the sky 
background, would roughly change as S/N} oc 1/nt. 
A way of comparing depths across different systems is 
the S/Nb in a fixed width band (obtained by com- 
bining all the individual filters included in that band): 
S/Nb oc 1/^/nJ. Therefore, by increasing the number 
of filters we would expect the effective limiting magni- 
tude at a fixed S/N level to diminish quite drastically, 
as mu m — const+2.5 log(^/nJ), e.g. equivalent to the 
loss of a full magnitude going from 4 to 25 filters. 

Fig. 2 describes how the 80% completeness magni- 
tude limit behaves for each of the filter systems. We 
see that for contiguous filters, the completeness depth 
sinks fast for rif < 8, and that the optimum number 
of filters is rif ~ 12, after which the effective complete- 
ness magnitude decreases, but much more slowly than 
expected from the change in the photometric limiting 
magnitude. This shows that for systems with low n/, the 
color-redshift degeneracies introduced by an insufficient 
wavelength resolution dominate over the improvement in 
S/N achieved by the increased filter width. 

Fig. 3 shows what happens when we add moderately 
deep near-IR observations with 5<r limiting (Vega) mag- 
nitudes of J = 22.4, H = 21.2, K = 20.4. There is a 
very significant, almost ~ 0.4 mag increase in the com- 
pleteness magnitude, and the behavior of the low n/ 
systems relatively improves, but still the most efficient 
overall performer is a logarithmically increasing band- 
width, half-band overlapping system with 9 filters, which 
reaches a completeness limit ~ 0.7 mag deeper than a 
typical 4 — 5 filter system with the same exposure time, 
while having a 10% better accuracy. 

Another obvious quantity to consider is the accuracy 
of the photometric redshifts for the high odds sample, 
estimated using the rms of the quantity Az/(1 + z) = 
[z— Zb)/(l + z), plotted in Fig. 4. Here we see that, as ex- 
pected, the redshift precision quickly and monotonously 
improves with rif and that adjacent filter systems per- 
form much better than overlapping ones. From Fig. 3 
we can see that a adjacent system with rif — 20 reaches 



Fig. 3. — Same as previous figure, but including moderately 
deep near-IR observations (see text for details) 
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Fig. 4. — Dependence of the rms of quantity (z — Zb)/(1 + z) 
for those galaxies with Odds > 0.99 as a function of the number 
of filters for the four types of filter system considered in the paper 
and including near-IR observations (see text for details). 
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a completeness depth similar to traditional systems with 
rif = 5, but an accuracy significantly better: 0.015(1 + z) 
vs 0.04(1 + z). 

4. THE ALHAMBRA SURVEY 

Considering the above results, plus additional require- 
ments on emission line detection, the ALHAMBRA (Ad- 
vanced, Large, Homogeneous Area, Medium Band Red- 
shift Astronomical) survey decided to use a constant- 
width, non overlapping filter system, complemented with 
near-IR observations. ALHAMBRA had an additional 
requirement, the detection of a large fraction of galax- 
ies with emission lines, which favored the choice of 310A 
filters. The ALHAMBRA 3cr rest-frame detection limits 
for a typical ABw 23 galaxy are EW(Ha)> 28A out to 
z « 0.45, and EW(OII)> 16A out to z w 1.55. From 
comparison with Hippelein et al. 2003, ALHAMBRA ex- 
pects to detect « 50% of the Ha emitters at z k 0.25, 
and w 80% of the Ollobjects to z « 1.2; since 80% of 
field galaxies at those redshifts have detectable emission 
lines (Tresse & Maddox 1998) we expect to detect lines 
for a large fraction of our whole sample. 
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The survey is imaging 4 sq. degrees with the camera 
LAIC A at the Calar Alto 3.5m telescope and also ob- 
taining deep JHK observations with Omega2000 at the 
same telescope. The survey characteristics, scientific goals 
and preliminary results are described in detail in Moles 
et al. (2008). A good test of the simulations presented 
in this paper is a comparison with ALHAMBRA. The 
mock catalogs predict that, with 20 filters, ALHAMBRA 
should be able to reach a precision of Sz/(1 + z) 0.014 
for / < 24 galaxies. Preliminary results show that the 
measured redshift error (Moles et al. 2008) is similar or 
less than 0.015, supporting the validity of the simulations 
presented in this paper. 

5. CONCLUSIONS 

We explore the performance of four different uniform 
filter systems with constant and logarithmically increas- 
ing (AA oc A) widths, and with half-width AA/2 over- 
laps or just a minimal overlap corresponding to the filter 
wings, and use, as a measure of survey effective depth, 
the fraction of galaxies with a compact, unimodal prob- 
ability redshift distributions as a function of magnitude. 
Our simulations employ a realistic input catalog, based 
on HDF photometric redshifts, and correspond to a com- 
bination of total exposure time and telescope area of 270 
hrs m 2 . We find that traditional 4 — 5 optical filter sys- 
tems clearly underperform, both in terms of complete- 
ness magnitude limit and precision, systems with rif > 8 
filters. 

Our results suggest that for low rif, the effect of color- 
redshift degeneracies dominates the advantages of in- 
creased photometric depth, and that even at higher rif, 
the effective photometric redshift depth decreases much 



more slowly with filter width than naively expected from 
the reduction in S/N. Adding near-IR observations in- 
creases the overall depth, alleviating color-redshift degen- 
eracies and improving the relative performance of low rif 
systems. However the optimum performance still comes 
from a system with 9 filters with logarithmically increas- 
ing bandwidth (constant resolution) and half-band over- 
lap, which reaches ~ 0.7mag deeper, with 10% better 
redshift precision, than 4 — 5 filter systems. For many 
scientific applications, which require both precision and 
depth, the use of > 15 medium band filters is clearly 
advantageous. A system with 20 constant-width, non- 
overlapping filters reaches only ~ O.lmag shallower than 
4 — 5 filter systems, but has a precision almost 3 times 
better, 5z = 0.014(1 + z) vs Sz = 0.042(1 + z), as a prac- 
tical implementation of such a system, the ALHAMBRA 
survey, shows. 

Since it is well known that color-redshift degeneracies 
worsen with magnitude depth, it can be expected that 
the relative decoupling between photometric depth and 
photometric redshift depth described here will be more 
significant for surveys which reach fainter limits than 
those considered in our simulations, and less important 
for shallower observations, where the color/redshift de- 
generacies are less of a problem. In any case, future 
projects will have to seek an optimum number of filters 
based on their particular observing parameters and sci- 
ence goals. 
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